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We consider a U(l)-invariant nonlinear Dirac equation, interacting with itself via the mean field mechanism. 
We analyze the long-time asymptotics of solutions and prove that, under certain generic assumptions, each 
finite charge solution converges as t — > ±00 to the two-dimensional set of all "nonlinear eigenfunctions" of the 
form <p(x)e~' tu ' t . This global attraction is caused by the nonlinear energy transfer from lower harmonics to the 
continuous spectrum and subsequent dispersive radiation. 

The research is inspired by Bohr's postulate on quantum transitions and Schrodinger's identification of the 
quantum stationary states to the nonlinear eigenfunctions of the coupled U(l) -invariant Maxwell-Schrodinger 
and Maxwell-Dirac equations. 



I. INTRODUCTION. BOHR'S TRANSITIONS AS GLOBAL ATTRACTION 



In the present paper we continue our research on the global attraction to solitary wave solutions in U(l)-invariant 
dispersive nonlinear systems. Our aim is a dynamical interpretation of "quantum jumps" postulated by N. Bohr in 1913 
to explain the stability of the Rutherford model of the atom (1911). 

The long time asymptotics for nonlinear field theory equations, and in particular the nonl inear K lein-Gord on equ ation, 
have been the subject of inten sive research, starting with the pioneering papers by Segal ISeg63ll . Strauss IStr68il , and 
Morawetz and Strauss lMS72ll . where the nonlinear scattering and the local attraction to zero solution were considered. 
The orbital st ability of solitary wave solutions of nonlinear Schrodinger and Klein-Gordon equations is well-understood 
(see llGSS87ln . and there is presently an active research on local attraction to solitary waves, or asymptotic stability 



I BS03LICuc08ll . 



Traditionally, the Dir ac equ ation has presented difficulties. The existence of solitary waves in the n onlinear Dirac 
equation was proved in ICV86ll . The existence of solitary waves in the Maxwell-Dirac syste m is pr oved in jEGS96ll . The 
stability of solitary waves with respect to a particular class of perturbations was analyzed in | SV86]. The spectral stability 
of small ampli tude solitary w aves of a particular nonlinear Dirac eq uation i n one dimension, known as the massive Gross- 
Neveu model IGN74 IlG75I1 . has been confirmed numerically in ICP06ll . Neither the orbital stability nor asymptotic 
stability of solitary waves is presently understood in the context of the Dirac equation. 

The next question is related to the structure of a global attractor of all finite energy solutions . Global attraction to static, 
station ary solutions in dispersive systems without U(l) symmetry was first established in lKom95l IKSK971 lKom99l 
iKSOOll . The attraction of any finite energy solution to the set of all soli tary waves in t he c ontext of the Klein-Gordon 
equation coupled to one and to several nonlinear oscillators was prov ed in jKK07llKK09ll . In llKK09ll . we generalized this 
result for the Klein-Gordon field coupled to several oscillators. In iKKOcl . a similar result was generalized to a higher 
dimensional setting for the Klein-Gordon equation with the mean field interaction. In the present paper, we consider the 
structure of the attractor for the Dirac equation, in the model with the mean field interaction. We will show that under 
rather general assumptions the attr actor co incides with the set of all solitary waves. 

According to Bohr's postulates ]Bohl3ll . an unperturbed electron runs forever along certain stationary orbit, which we 
denote \E) and call quantum stationary state. Once in such a state, the electron has a fixed value of energy E, with the 
energy not being lost via emitted radiation. Under a perturbation, the electron can jump from one quantum stationary 
state to another, \E-) 1 — > \E+), radiating (or absorbing) a quantum of light. Bohr's postulate suggests the dynamical 
interpretation of Bohr's transitions as a long-time asymptotics 

ip(t) — >\E±), t^±oo (1.1) 
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for any trajectory tp(t) of the corresponding dynamical system, where the limiting states \E±) depend on the trajectory. 
Then the quantum stationary states should be viewed as points of the global attractor. 

At first glance, the global attraction (II. It seems incompatible with the energy conservation and time reversibility of 
Hamiltonian systems. We intend to verify that such asymptotics in principle are possible for nonlinear Hamiltonian field 
equations. In this paper, we verify this asymptotics for equations of Dirac type. 

Developing de Broglie's ideas, Schrodinger identified quantum stationary states of energy E with the solutions of type 

ip(x,t) = e- lEt/7l ^(x), i€R 3 . (1.2) 

Then the Schrodinger equation 

h 2 

iW t ip = Hip := -—Aip + V(x)ip, V = ip(x,t) € C, x e K 3 , (1.3) 
2m 

becomes the eigenvalue problem Eip(x) = Htp(x). One of the original Schrodinger's ideas | Sch26 a|] was to identify the 
integers in the Debye-Sommerfeld- Wilson quantum rules with the integers arising in the eigenvalue problems for PDEs. 

For the case of the free particles, this identification agrees with the de Broglies wave function ip(x, t) = ^P-xfhg-tEtfh^ 
where p £ M 3 is the momentum of the particle. For the bound particles in an external potential, the identification (II.2b 
reflects the fact that the space is "twisted" by the external field, while the time remains "free". 

Thus, the reason for Schrodinger's choice of quantum stationary states in the form dl.2| ) seems to be rather algebraic. 
At the same time, the attraction ( II. U suggests the dynamical interpretation of the quantum stationary states as asymptotic 
states in the long time limit, that is, 

ip(x, t) ~ cf> u± (x)e- iw±t , t -> ±oo. (1.4) 

Such asymptotics should hold for each finite charge solution. The asymptotics of type dl.4t are generally impossible for 
the linear autonomous equation ( II. 3b because of the superposition principle. The asymptotics would mean that the solitary 
waves il.2\ form the global attractor for the coupled nonlinear Maxwell-Schrodinger and Maxwell-Dirac systems. 

An adequate description of "quantum jumps" in an atom requires that we consider the equation for the electron wave 
function (Schrodinger or Dirac equation) coupled to the Maxwell system which governs the time evolution of the four- 
potential A^{x, t) = {jp {x, t), A(x,t)). The corresponding Maxwell -Schrodi nger system was initially introduced by 
Schrodinger in ISch26bll . Its global well-posedness was cons idered in The results on local existence of solu- 

tions to the Dirac-Maxwell system were obtained in ]Bou96ll . The coupled systems are U(l)-invariant with respect the 
(global) gauge group Up(x), A fJ -(x)) H> (<p(x)e~ 10 , A^(x)), 9 € K. Respectively, one might expect the following natural 
generalization of asymptotics ( II-4b for solutions to the coupled Maxwell-Schrodinger or Maxwell-Dirac equations: 

{^(x,t),Ai*(x,t)) ~ ((p ± (x)e- iw±t ,A^(x)), i^±oo. (1.5) 

The asymptotics dl-5b would mean that the set of all solitary waves forms a global attractor for the coupled system. The 
asymptotics of this type are not available yet in the context of coupled systems. 



II. MODEL AND THE RESULT 



In the present paper we consider the Dirac equation with the mean field interaction: 



iip(x,t) = (-i^atjdj +(3m)<ip + p{x)F((p,ip(;t))), i/»eC 4 , x € R 3 , n > 1, 
j'=i 



(II. 1) 



where the Hermitian 4x4 matrices otj and (3 satisfy {otj, = 25jkl, {otj,/3} = 0, /3 2 = I. We write (p, tp(-, t)) 
/ p' (x)ip(x, t) d 3 x, where p is a spinor-valued coupling function from the Schwartz class: 

JR3 



p = 



e^(R 3 ,C 4 ), P ^o, p t = [p 1; ...,p 4 



Assumption II.l. F(z) = — VRcz.imzU(z) with U(z) — J2 Uk\z\ 2k , where <ER,p> 2, and u p > 0. 

k=l 
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Since one has F(e lti z) = e F(z) for 9 £ E, z £ C, equation ( III. II ) is U(l)-invariant, where U(l) stands for the 
unitary group e l °, 9 £ R mod 27r. By the Nother theorem, the charge functional Q(tp) — J ffi3 \ip\ 2 d 3 x is conserved. 

Equation ( III. It is the simplest model sharing the following features with the Maxwell-Dirac system, which we consider 
the key for global attraction to solitary waves: The model has U(l)-symmetry; It is dispersive; It is nonlinear. The first 
feature allows for solitary waves, while the last two features are responsible for the "friction by dispersion" mechanism: 
The nonlinearity moves the perturbations into the continuous spectrum of the linearized equation, and thereafter the 
perturbations are dispersed to infinity. Let us also mention that the modelling of th e interac tion of the matter with gauge 
fields in terms of local self-interaction takes its origin back to at least as early as ISch5lll . where the Lorentz-invariant 
nonlinear Klein-Gordon equation originally appeared. 

Definition II.2. The space of states of finite charge is 56 = L 2 (R 3 , C 4 ), with the standard norm denoted by || • ||^2. 

Equation ( III. Il l formally can be written as a Hamiltonian system with the phase space 56 ', ip(t) = J DH(ip), where 
DW is the variational derivative with respect to Re Im tp^, [i = 1, . . . , 4, of the Hamiltonian 

^(V0 = 7r/ i>H-ia j d j +0m)^^x + U({p,il;)), = $ 4 ]- 

^ ./m?3 



Let = OL^j + (3m, £ £ R 3 . Since is self-adjoint and £^(£) 2 = £ 2 + m 2 , its eigenvalues are ±-\/£ 2 + m 2 

Define the projectors II ± (£) = i ( 1 ± ) onto the corresponding eigenspaces, and denote p± (£) := IT± (£)p(£). 

V v ? 2 + m2 ' 

Assumption II.3. (/) For any A > 0, /3+(£) and /3_(£) do not vanish identically on the sphere |£| = A. 
(ii) The function 



i \ I pHt)(u + * j Z j +l3m)p(Z) d 3 f 

^ = ' ^-e-m 2 ( } 



is nonzero for all a; £ [— m, m], except perhaps at one point which will be denoted uj a £ [—m,m]. 

Remark II.4. Note that limits <r(m — 0), <r(—m + 0) exist. 

Let us give examples of such functions p £ ^(R 3 , C 4 ). We take a spherically symmetric function p such that />(£) 7^ 
for £ 7^ 0, with /3(£) being an eigenvector of /3 (since /3 is Hermitian and (3 2 — 1, its eigenvalues are ±1). Then 

\P±(0\ 2 = |pt(C)(l ± ^f±5)p(6> hence 

f \p±(t)\ 2 d 2 n> [ p(o t (i-- 7 ^=)p(0d a n>o ) a>o. 

J|£|=A •'IfNA + m 

The integral of the term with cty£j dropped out since p{6) is spherically symmetric. 

Let us now consider cr(u>). Due to p being spherically symmetric, ct^-term cancels out from the integration in ( III.2I >. 
Since we require that be an eigenvector of /3, one has p{(6)* (lo + (3m) p{6) = (w±m)/5(£)*/5(£), hence the expression 
under the integral in ( |II.2t is sign-definite for all lj £ [—m, m] except at lj = m or at lj = — m. Therefore, (t(lj) does not 
vanish on [—m, m] except at one of the endpoints. 

Definition II.5. The solitary waves of equation (III. lb are solutions of the form ip(x,t) = 4> LJ (x)e~' lut , where u £ R, 
<fiu;(x) £ 56 ' . The solitary manifold is the set S = {(f)^: 4> UJ {x)e~ lu ' t is a solitary wave} C 56 . 

Denote C + = {w e C: Imw > 0}. For u £ C + U (— m, m), introduce 

m U ) = + E(x, U ) = [£(-, a,)] • (H.3) 

lj z — — m z 5 

E(-, lj) is an analytic function of lj £ C + with the values in ,5^(R 3 , C 4 ). One can show that for any x £ R 3 , E(x, lj) can 
be extended to the real line lj £ R as the boundary trace (in the sense of tempered distributions): 

E(x,lj) = lim E(x,u + ie), uel. (II.4) 

e— >0+ 

One has E(-, lj) £ 56 for lj £ (— m, m). For lj = ±m, I7(-, cj) € 56 only when p±(0) = 0. 
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Proposition II.6. Let Assumptions \n.l\ and\ lL3\ be satisfied. Then there are no nonzero solitary wave solutions to ( III. It 
for uj (f. [—to, to]. For uj £ (—to, m), there is a nonzero solitary wave (j) UJ (x)e~ , ' ult if and only if there is C £ C\0 such 
that 

-a{w)U'{Ca(w)) = l, (11.5) 

with U(-) from Assumption UI.l] One has (f>uj(x) — CS(x,oj). For uj — ±to, </> w (x) £ 3£ only when /5±(0) = 0. 

Remark II.7. Due to Assumption lII.il the set S from Definition III. 51 is invariant under multiplication by e tB , 6 £ R. 
Generically, the solitary manifold S is two-dimensional. It may contain disconnected components. 

Remark II.8. Since F(0) = 0, the zero solitary wave is an element of S, formally corresponding to any 

Remark II.9. It is possible that for some uj equation dll. 5b has nonzero roots Cj £ C of different magnitude; then there 
are solitary wave solutions corresponding to each of these roots. 

Definition 11.10. Fix e > and \ € C^ mp (R 3 ), < x < 1> x( x ) = 1 for M < 1> with suppx in the ball of radius 
2. Denote by W the space X endowed with the metric Hulls' := X>.ReN 2~ R \\x( x / fyi-'i^Wn-^ where the norm in the 
Sobolev space H s is given by || -011 h> = \\( m2 — A) s / 2 -0|| i 2. 

The Sobolev Embedding Theorem implies that the embedding S£ C & is compact for any e > 0. 

Theorem 11.11 (Main Theorem). Assume that the nonlinearity F(z) is as in Assumption 1/7.71 and that the coupling 
function p(x) satisfies Assumption 177.51 Then for any ipo £ SC there is a global solution ip(t) £ C(R, to equation 
( III. II ), such that V'| t=0 — V'o. which converges to S in the space *ty : 

lim dist ^ (^(i), S) = 0, where dist ^(V', S) := inf II?/; — sll^. 
t— >±oo ses 




FIG. 1 : For t — > ±oo, a finite charge solution ip(t) approaches the global attractor which coincides with the set of all solitary waves S. 

We work in the L 2 setting because of the absence of existence results for finite energy solutions. As a result, the conver- 
gence to the attractor holds just below L 2 , in local H _E -norm. The convergence to S is due to the "friction by dispersion": 
the nonlinearity carries the excess energy into the continuous spectrum, where it is subsequently dispersed to infinity, and 
the remaining part of the s olution settles to one of the nonlinear eigenstates (or zero). By the Titchmarsh Convolution 
Theorem (see llTit26ll and lYos80|P . this process can only stop when the time-spectrum of the solution consists of the 
single frequency, that is, when the limiting solution coincides with a nonlinear solitary wave (a nonlinear Schrodinger 
eigenstate). It is only for such states that the energy transfer from lower harmonics into the continuous spectrum (the 
spectral inflation) is absent. 

Let us comment on our methods. We follow the path developed in llKK07[ IkK081 IkK 091 : we prove the absolute con- 
tinuity of the spectral density for large frequencies, use the compactness argument to extract the omega-limit trajectories, 
and show that these trajectories have the finite time spectrum. Then we use the Titchmarsh Convolution Theorem to 
pinpoint the spectrum to a single frequency. Let us note that the Titchmarsh theorem is only applicable if we assume that 
F admits a real-valued polynomial potential (Assumption lll.ll i. 



III. CONCLUSION 



We demonstrated that in a particular nonlinear dispersive equation, whose linearization is the Dirac equation, the conver- 
gence to Schrodinger states for long positive and negative times takes place as a consequence of "friction by dispersion", 
so that the solution settles to one of the nonlinear eigenstates (or to zero). We conjecture that a similar consideration will 
allow to prove convergence to nonlinear eigenstates in systems such as coupled Maxwell-Dirac equations, providing a 
dynamical description of Bohr's "quantum jumps". 
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